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Abstract 

The Next -to-Minimal Supersymmetric Standard Model (NMSSM) includes a Higgs 
iso-singlet superfield in addition to the two Higgs doublet superfields of the minimal 
extension. If the Higgs fields remain weakly coupled up to the GUT scale, as 
naturally motivated by the concept of supersymmetry, the mixing between singlet 
and doublet fields is small and can be treated perturbatively. The mass spectrum 
and mixing matrix of the neutralino sector can be analyzed analytically and the 
structure of this 5-state system is under good theoretical control. We also determine 
decay modes and production channels in sfermion cascade decays to these particles 
at the LHC and pair production in e + e~ colliders. 



1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) p!J |2j opens the path to the 
analysis of supersymmetric theories. Arguments have been advanced however that suggest 
extensions beyond this minimal version. One well-motivated example is the Next-to- 
Minimal Supersymmetric Standard Model (NMSSM) |3] in which an iso-singlet Higgs 
superfield S is introduced in addition to the two iso-doublet Higgs fields H Ut d incorporated 
in the MSSM to generate electroweak symmetry breaking. Such an extension offers a 
possible solution of the \i problem, generating in a natural way, a value of the order of 
the electroweak breaking scale v; this is achieved by identifying /i, apart from the 0(1) 
coupling, with the vacuum expectation value of the scalar component S of the new iso- 
singlet field. [For a recent summary of this construct see Ref.[4j; a useful code has been 
made available in Ref.|3].] 

The superpotential of the NMSSM includes, besides the usual MSSM Wy Yukawa 
components, an additional term, which couples the iso-singlet to the two iso-doublet 
Higgs fields, plus the self-coupling of the iso-singlet: 

W = Wy + XS(H u H d ) + ^kS 3 (1) 

The two parameters A and k are dimensionless. By demanding the Higgs fields remain 
weakly interacting up to the GUT scale, the two couplings are bounded at the electroweak 
scale by the inequalities A, k < 0.7. While the scalar Higgs sector includes several soft 
supersymmetry breaking parameters, the Lagrangian of the gaugino/higgsino sector is 
complemented only by the familiar SU(2) and U(l) gaugino mass terms. As a result, the 
parameter space of the neutralino sector is much less complex than the Higgs space. 

The superpotential without the singlet self-coupling, i.e. n = 0, incorporates a Peccei- 
Quinn (PQ) symmetry: {H U {1), H d (l), S{-2), Q(-l), U(0), D{0\, E(0)}. Q 

and L are the quark and lepton SU(2) doublet superfields, while U, D and E are the 
up- and down-quark and lepton SU(2) singlet superfields respectively. The integer of 
each parenthesis indicates the PQ charge of the corresponding superfield. The sponta- 
neous breaking of this symmetry by the non-zero vacuum expectation value v s / V2 of the 
scalar S field gives rise to a massless Goldstone boson. However, when k ^ 0, the mass 
is lifted to a non-zero value by the self-interaction of the 5* field. Still, a discrete Z 3 
symmetry is left which would lead to the formation of domain walls in the early Universe. 
This problem can be tamed by introducing new interactions of the inverse Planck size 
that, however, do not affect the low-energy effective NMSSM theory jH]. 

In contrast to the Higgs sector, masses and mixings in the chargino system are not 
affected by the singlet extension. [Of course new decays such as S — > x + X~ or X + X§H + 
may be possible if allowed kinematically.] 

So far the supersymmetric particle spectrum of the NMSSM has received only little 
attention in the NMSSM literature, Refs.jH], 0-0 ■ In this report we attempt a system- 
atic analytical analysis of the neutralino system. In contrast to the MSSM where exact 
solutions of the mass spectrum and mixing parameters can be constructed mathemati- 
cally in closed form, this is not possible any more for the NMSSM in which the eigenvalue 
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equation is of 5th order, not allowing closed solutions. However, since the coupling be- 
tween singlet and doublet fields is weak, Af /a/2 ~ O(10 2 ) GeV, compared with the typical 
super symmetry scale Mi j2 and /i = \v s /\/2 ~ O(10 3 ) GeV, a perturbative expansion of 
the solution gives rise to a good approximation of the mass spectrum while the magnitude 
of the matrix elements in the mixing matrix is at least qualitatively well understood. The 
usefulness of a perturbative expansion has also been noticed in Ref. [OJ; however, here, 
extending the Higgs analysis in Ref. [4., we work out this approach systematically for all 
facets of the NMSSM. 

While k plays a crucial role in the Higgs sector, it is less crucial for the neutralino 
system. The size of kv s , with v s < 15 v to maintain a link with the electroweak scale, just 
determines the singlino mass before modified by mixing effects. Once masses and mixings 
are determined, the couplings of the neutralinos to the electroweak gauge bosons and to 
scalar /fermionic matter particles are fixed. Decay widths and production rates of the five 
neutralinos can subsequently be predicted for squark cascades at the LHC [TU] and e + e~ 
annihilation at prospective linear colliders [TT] . 

The report is organized as follows. In Sect. 121 we describe the neutralino sector of 
super symmetric models in which the pair of Higgs doublet superfields is augmented by an 
additional iso-singlet field. In Sect.|H]we show how, for a naturally expected weak coupling, 
the properties of the four standard neutralinos are modified; moreover the properties of 
the fifth neutralino, the new singlino-dominated state, are calculated. All these spectra 
and mixings are pre-determined analytically before the surprisingly good quality of the 
weak-coupling expansion is demonstrated by comparison with numerical solutions. In 
this way we achieve a satisfactory theoretical understanding of the system. In the limit 
of large gaugino mass parameters M± t 2 compared with the higgsino mass parameter /i, or 
vice versa, the MSSM part can be easily diagonalized analytically and a clear and simple 
picture of the entire system emerges. The section is concluded by a lovely toy model 
in which we set Mi = M 2 and tan/3 = 1; this set allows us to solve the system exactly, 
leading to transparent closed expressions for the neutralino mass spectrum and the mixing 
parameters. A sample of decay widths and production cross sections for the neutralinos 
is presented in Sect. EJ The results are summarized in Sect. 03 and technical details of 
the diagonalization procedure for the 5x5 neutralino mass matrix are described in the 
Appendix. 

2 The NMSSM Neutralino Sector 

2.1 The NMSSM neutralino mass and mixing matrix 

The Lagrangian of the neutralino system can be derived from the superpotential defined 
in Eq.(£Q), complemented by the SU(2) and U(l) mass terms in the soft supersymmetry 
breaking Lagrangian. After breaking the [electroweak] symmetry spontaneously by intro- 
ducing non-zero vacuum expectation values of the iso-doublet and singlet Higgs fields, 

W = -j= sin/3 ( ° ) , (H d ) = -L cos/3 ^ ) , (S) = v s /V2 (2) 
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the Higgs-higgsino mass parameter 



(3) 



is generated and, subsequently, the neutralino mass matrix 

(MX 

with a hierarchical structure as analyzed in the Appendix, can be written in detail as: 
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This 5x5 mass matrix is constructed from the standard 4x4 MSSM neutralino mass 
matrix M. in the upper left corner, the mass term of the higgsino component S of the 
singlet superfield S, 



2kv< 



and the mixing between doublets and singlet parameterized by 



/'a 



Xv 



(5) 



(6) 



The mass matrix J\4 5 is defined in the group basis (B, W 3 , H%, H®, S). As usual, Mi and 
M 2 are the soft SUSY breaking U(l) and SU(2) gaugino mass parameters, tan/3 is the 
ratio of the vacuum expectation values of the two neutral SU(2) Higgs doublet fields (as 
defined in Eq. (J2J)), sp = sin (3, cp = cos (3, and s\y, cw, tw are the sine, cosine and tangent 
of the electroweak mixing angle 8w 

Since the neutralino mass matrix (@J) is symmetric and real, it can be diagonalized by 
an orthogonal matrix V 5 . The mass eigenvalues are real but not necessarily positive. They 
can be mapped onto positive values by supplementing the rotation matrix to iV 5 = $ 5 U 5 
with the diagonal phase matrix {& 5 )ki = $ki in case of positive (negative) eigenvalues 
so that iV 5 * A^5 N 5 ^ is positive diagonal. The physical neutralino states Xi [i = 1 — 4] are 
ordered according to ascending mass values while x° is the predominantly singlino state.* 
They are mixtures 



.vj(/i.U' :i . 



Xi 



Hd, H u , S) 3 



1-5] 



(7) 



*Note that the ordering of the masses according to ascending values is accomplished easily after the 
diagonalization process is finalized. For the intermediate steps it is however convenient to use the indices 
i = 1,2,3,4 for the former MSSM type states and i = 5 for the additional state originating from the 
singlino field as suggested by the structure of M5 in Eq.(4). 
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of the U(l), SU(2) gauginos, the doublet higgsinos and the singlino. 

The unitary matrix N 5 defines the couplings of the mass eigenstates Xi to other par- 
ticles. For the neutralino production processes it is sufficient to consider the neutralino- 
neutralino-Z vertices 



(x°l\Z\x jL ) 

(x° R \z\x° jR ) 



2c w 
9 



(N%N* 3 - n?:n* a ) 



2cw 



and the fermion-sfermion-neutralino vertices 



(x° R \h\h) 

C& I /*!/*> 



-V2— 
c w 



2gQf s w N, t 



llN?*c w + {Q f - Il)N* 8w 



a 



(9) 



The coupling g is the SU(2) gauge coupling, // is the SU(2) isospin 3-component and Qf 
is the electric charge of the fermion /. In Eq. (JHJ) the coupling to the higgsino component, 
which is proportional to the fermion mass, has been neglected for "light flavors". The 
more involved Higgs couplings to the neutralinos are listed in detail in Sect. EJ 



2.2 NMSSM parameter range 

In contrast to the Higgs sector only two additional parameters A and k are introduced in 
the NMSSM neutralino sector as compared to that of the MSSM including fi. Assuming 
that the fields remain weakly interacting up to the GUT scale, the two couplings are 
bounded at the electroweak scale by the inequality 

VA 2 + ^ < 0.7 (10) 

Moreover, the renormalization group (RG) evolution of the couplings points to k < A 
as preferential target domain if the evolution starts from a random distribution of the 
couplings kjj, A(7 < 2n at the GUT scale 0j. 

While v = \fv^ + v\ = 246 GeV is fixed by the Fermi coupling Gp, the parameter v s 
should be expected in the same range, 

v s <15v (11) 

in compliance with the arguments for introducing the NMSSM. A RG analysis of the entire 
set of parameters shows that a low value of tan/3 is favored 0j. Current experimental 
analyses of tan/5 assume MSSM relations for the couplings; they are modified in the 
NMSSM and the results in this extended scenario are less restrictive. 

Since the size of the doublet -singlet mixing is set by = Xv/y/2, the mixing in- 
teraction is expected to be small^ compared with the standard supersymmetry scales, 

* [i\ is expected to have a lower limit from cosmological arguments; private communication with U. 
Ellwanger, see also Ref.|12l. For too small a value of i.e. very much below the typical scale v/y/2, 
the amount of cold dark matter may exceed the measured value of Ocdm ~ 0.25; detailed analyses are 
not available yet. 
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[i = \v s /\/2 and/or M ii2 for which values < 0(1 TeV) are anticipated. As a result, 
transparent expressions can be found by performing a systematic expansion for small 
mixing between the gauginos/doublet higgsinos and the singlino, measured by the small 
size of the parameter /i^ relative to the other parameters in the mass matrix. 

In summary, at tree-level the NMSSM neutralino sector described above has six 
free parameters which we choose as fi K and fi\ in addition to the MSSM parameters: 
{{Mi, M 2 , tan/?, /x}; Sometimes it is convenient to re-express fi, fi\ and \i K in 

terms of A, k and v s . The spectrum of the NMSSM neutralino sector will now be ana- 
lyzed in detail. 



3 NMSSM Small-Mixing Scenarios 

In general the diagonalization of the 5x5 NMSSM mass matrix A4§ cannot be performed 
analytically in closed form. However, if the doublet-singlet coupling is weak, an approxi- 
mate analytical solution can be found after the 4x4 MSSM submatrix M. is analytically 
diagonalized following the elaborate standard procedures in Ref . [13J . 

The orthogonal matrix V 5 which transforms M.§ to the diagonal mass matrix 
M.® = diag [mi, . . . ,m^\ is conveniently split into a matrix V diagonalizing the 4x4 
submatrix A4 and a matrix performing subsequently the block diagonalization of the 
4x4 and lxl submatrices. After the block-diagonalization, the upper left MSSM mass 
matrix Ai D = diag [m^mj, m^, rh^} needs not be re-diagonalized for small doublet-singlet 
mixing, as proved in the Appendix. The final result for the orthogonal matrix V 5 may be 
written in the simple form: 

§(VT)(Vr) T (VT) \ ( V 



4x4 



i - Uvv) T (vr) 



o li 



xl 



(12) 



The doublet-singlet 4-component mixing vector T can be expressed in terms of the gaug- 
ino/higgsino parameters as 

/ M^fimz s w c 2 /3 \ 

-M[ fj, m z c w c 2f 3 
M[ M' 2 {ii c p - fi K sp) - M[ 2 m 2 z s p 
\ M[ M' 2 {ii sp - fi K cp) - M[ 2 m 2 z cp ) 



det{M -fi K ) 



(13) 



with the abbreviations 

M[ = Mi - fi K , M', .\/, - fji K , MS = M[' } c 2 w + M 
and the determinant 

det(M - aO = M[M' 2 (iil - ii 2 ) + M[ 2 {fxs 2 p + fi K ) m 



The mixing with the singlet alters the MSSM mass eigenvalues rrij 

r2U 



0(e 2 y, and correspondingly the singlet mass 



m 5 = ji K 



(14) 

(15) 
4] to 

(16) 



*Note however that small mass differences 



TO5I may enhance the mixing effects. 
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The shifts are given as 



m i = m l + - 1 . {VXft [i = l-4] 
mj — m 5 

m 5 = rh 5 - (VX)l (17) 

/ — J rrii — m 5 

with the 4-component vector X = fj,\(0, 0, — sp, —cp). [The eigenvalues are not necessarily 
ordered sequentially, and, if some of them are negative, the additional phase rotation 
transforms them to positive physical masses.] Even for small mixing, the 5th eigenvalue 
m 5 may differ significantly from the singlino mass parameter m 5 = fi K if k is small. 
However, even though the relative shift may be large, the absolute shift remains small, of 
second order. Trivially, the eigenvalues fulfill the spur formula 

5 

J2^ = M 1 + M 2 + fi K (18) 
i=i 

which is independent of the parameters /i and 

The doublet-singlet mixing generates a singlino component in the wave functions of 
the original MSSM neutralinos x° [i — 1, • • • , 4] of the size 

4 
3=1 

linear in the mixing parameter to first approximation as expected for off-diagonal ele- 
ments. Reciprocally, the singlino component in the wave function of is reduced to 

^i-^Er? (20) 



i=i 



differing from unity only to second order in the mixing as expected for diagonal elements. 

As long as the mixing parameter is significantly smaller than the other parameters, 
we find that the approximation works remarkably well, as demonstrated in Fig. ^ As 
an example both the exact numerical solutions and the approximate solutions for the 
neutralino masses are shown as a function of fi\ for a favored parameter set P of broken 
PQ symmetry, //„ = 120 GeV with M 1 = 250 GeV, M 2 = 500 GeV, fi = 170 GeV and 
tan/? = 3. The exact and approximate solutions agree rather well as long as fi\ is less 
than about 80 GeV, as the mixing corrections are of second order in 

In Fig.|2]the exact numerical solution (solid) and the approximate solution (dashed) are 
compared for the gaugino/higgsino and singlino components, {|(A r5 ) 5 i|, |(A^ 5 ) 53 |, |(-/V 5 ) 55 |}, 
of the lightest singlino-dominated neutralino as a function of fi\ for the same parameter set 
P. Since the matrix V 5 is in general linear in the mixing term fi\, the approximate solution 
differs from the exact solution already for smaller values of fi\ in the reference point P 
in which m 5 = /i K is quite close to the higgsino parameter /i, though the characteristic 
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Figure 1: The exact numerical solution (solid) and the approximate solution (dashed) 
for the masses of the five neutralino states in the NMSSM as a function of fix for the 
parameter set P = = 120 GeV, M 1 = 250 GeV, M 2 = 500 GeV, fi — 170 GeV, 
tan/? = 3}. The ordering of the mass spectrum is m 5 ,mi,m 2 ,m 3 , and m 4 in increasing 
mass, i.e. the state x° is the lightest neutralino for the given parameter set. 



features remain valid up to /i\ ~ 40 GeV. 

To fully exhaust the potential of our analytical method we perform the complete 
NMSSM diagonalization for the two standard limits analyzed in general within the MSSM: 
3> |/i| and vice versa, both complemented of course by small doublet-singlet mixing 
fix < max{Mi >2 , 

3.1 Small singlino mass parameter 

The first special analysis should be performed for small singlino mass parameter fi K , which 
implies a slightly broken PQ symmetry kC 1 as favored by the RG flow of this coupling 
in grand unified theories. Due to the small doublet-singlet mixing the structure of the 
original MSSM neutralinos Xi [* = 1 — 4] is changed little while the properties of the 
5th neutralino Xs? the lightest for small fi K , are determined jointly by both the singlino 
parameter /i K and the mixing parameter 

3.1.1 Large gaugino mass parameters 

As a first example, we consider the case with large gaugino mass parameters, i.e. 
M li2 > |/i| > m z , fix- 
To begin, the 4x4 diagonalization matrix V defined in Eq. (|12j) can be parameterized 
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Figure 2: The exact numerical solution (solid) and the approximate solution (dashed) 
for the gaugino/higgsino and singlino components, {1-ZVfjJ, \N$ 3 \, | iVf 5 1} ; of the lightest 
singlino- dominant neutralino as a function of fi\ for the same parameter set P as in 
Figm 



up to second order according to standard MSSM procedure, cf. Ref . [T3] . as 

/V G \ /l ax , V x Wl axa \ 



V H J \ ~Vx ^2X2 J V itV/4 

The 2x2 7r/4 rotation = (1 — ia y )/V2 shifts the [34] off-diagonal elements [—fi, —fi] 
onto the diagonal axis [p,, —fi]. The second matrix, Vx, 

( -c p s w m z lM x s p s w m z /Mi \ 
V x = (22) 
\ c^c^mz/Mj -s f3 c w m z /M 2 J 

removes the mixing between the blocks of the two gaugino and the two higgsino states. 
The components Vq and Vjj diagonalize the gaugino and higgsino blocks themselves: 



2x2 



^ml/Mf 

4m|/M| 



_ 1 / (1 + M'UmlllMlMl \ 

H ~ 2x2 2 ^ Q (1 _ a ^ ^ m | /2M 2 M 2 J I J 

with Mf| = M\c\j + M|s^, respectively. Vg and V# relate to a diagonal form of the 
gaugino-higgsino mass matrix for large M 12 and fi. Their off-diagonal matrix elements 
are of second order and can be omitted consistently as they would effect the eigenvalues 
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only to fourth order. 

After these steps are performed, the 4x4 mass submatrix is diagonal and the complete 
symmetric mass matrix M5 takes the form 

/mi \ 

m 2 

m 3 +^\c- 



M 5 -> Mi 



m 4 



(24) 



\ +/i A C_ ~/iAC+ //« / 



where, in an obvious notation, zero elements are suppressed for easier reading, and 
c± = (c/3 ± sp)/y/2 is used as abbreviation. The MSSM neutralino mass eigenvalues are 
given by 



mi 
m 2 
m 3 = n 



m\ 
Mi 

™l 2 



M- 



12 



m4 = —ji 



2M 1 M 2 

Mia 



-m|(l - S2/3) 



(25) 



2MiM 2 

It remains to diagonalize .M5 by choosing the proper form of in V 5 . 

In the limit of large gaugino mass parameters, the doublet-singlet 4-component mixing 
vector T reduces to a simple expression 



(J* 



(0, 0, C/3, s p ) J 



(26) 



and the entire matrix V 5 can be written, up to second order, in the form 



/ 1 



V 5 



1 



^'1 + 820) 



4 M - 
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MA 







— — c_ 



— — C_i_ 



1 - iS- i 




(27) 



with zero's suppressed in the upper 4x4 matrix, and antisymmetric in the off-diagonal 
elements. 

The rotations lead eventually to the diagonal mass matrix, of which the mass eigen- 
values to the desired order are given by 



mi 



my 9 



m 



2 

Z J2 



m2 K M 2 + ^ 

m 5 ~ fi K + —s 2 p (28) 
/i 

[recall M\i = MiCyy + M 2 s'^ v }. For the ordering of the eigenvalues and the flipping of the 
signs to positive physical masses the previous general remarks apply. 

Two points should be emphasized explicitly. While the large gaugino masses mi, 7712 
are not affected by the singlino, it does affect the higgsino states 3,4 to second order. The 
singlino mass is also affected to second order; however the mixing term can be leading if 
the singlino mass parameter \i K is small. 

The mixing in the wave-functions is described by the components of T itself [since the 
4x4 matrix V deviates from unity only to second order in the small parameters of the 
order of the SUSY scales] : 

^(C/3 + S/3) 





A 4 


(°' °' ~^ C/3 ~ S ^ 




A 4 


(0, 0, Cp, 8p\ 


v 5 \ « 


f - 


A*A 

2^ 2 



(29) 

3.1.2 Large higgsino mass parameter 

As a second example, we consider the case with large higgsino mass parameter, i.e. 
\fi\ ^> M 1;2 ^ m z ,^\- This example is complementary to the previous case. 

The overall diagonalization 4x4 matrix V can be parameterized in the same form as 
that in Eq.(|2ip. The 2x2 matrix Vx describing the mixing between the two ensembles 
of the gaugino states and the higgsino states reads 



A« V ~ Cw S P ° w C P 



(30) 



leading to a block-diagonal mass matrix composed of a 2 x 2 matrix, depending on Mi 
and M2 with small corrections of the order of m 2 z / /1, and a 2 x 2 mass matrix, depending 
only on the higgsino parameter /1. The 2x2 blocks Vq and Vh in the gaugino and higgsino 
sector may be written 

V G w 1 2X 2 - - 
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|(1 - s 2 ^)m|//i 2 



(31) 



respectively, after the higgsino submatrix has been diagonalized by the standard R^/a 
rotation. 

These transformations diagonalize the 4x4 submatrix within the block-diagonal ma- 
trix M.§ of the same form as (|24j). of which the first four diagonal elements are given 
by 



ni 



fhx = Mi - —s 2 w s 2/ 3 



m 



m 2 = M 2 - —(? w s 2/ 3 
/i 



m. 



m 3 = jj + — (1 + s 2 p) 



m. 



(32) 



The mixing between the doublet-higgsino and singlino states is then described in an 
analytic form by a 4-component column vector 



r«-(o, 0, Cp , Sf3 f 

Z 1 



(33) 



mixing the singlino both with the gauginos and with the doublet -higgsinos. The entire 
matrix V 5 can be written up to second order in the form, with antisymmetric off-diagonal 
elements, 
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(34) 



and with the same abbreviations as before in Eq. (j24|) . 

The rotations lead eventually to a diagonal mass matrix, consisting of the mass eigen- 
values: 



mi » Mi - — s 2 w s 2 p 



rn 



m 2 « M 2 - — s 2/3 
m 8 «A*+g(l + « V ) + ^(l— V ) 



A* a , 



11 



m 4 « -fi - -^"l 1 _ s 2/3j 
Ma 




(1 + S 2j 9) 



(35) 



with apparent reciprocity in the MSSM subsystem between gaugino and higgsino param- 
eters in comparison with the previous case, but universal modifications from the doublet- 
singlet mixing. 

Correspondingly, to leading order the coefficients of V 5 involving the singlino index 5 
coincide with the elements of the doublet-singlet mixing matrix in Eq. (j29|) . 

3.2 Large singlino mass parameter 

In the alternative extreme, the PQ symmetry is strongly broken if k is large and, equiva- 
lently, the singlino mass parameter is large, i.e. fi K 3> fix, fi, M\p. This limit is not favored 
by the renormalization group flow from the GUT scale down to the electroweak scale but 
cannot be ruled out a priori on general grounds. The new fifth eigenstate, predominantly 
composed of the singlino, would in general be the heaviest state, mixed only weakly with 
the iso-doublets and, as a result, coupling weakly to electroweak gauge bosons and matter 
fields. 

Applying the approximation method described in the Appendix and the general in- 
troduction to this section, the neutralino mass matrix can be transformed into the 4x4 
and lxl block-diagonal form by inserting the mixing column vector 



in the V 5 matrix Eq.()12p. Note that the mixing column vector (|3*Bj) is directly proportional 
to the 4-component off-diagonal column vector of the mass matrix (jlj) unlike the column 
vector (|26|) for a small singlino mass parameter. 

From the general analysis it is apparent that the first four neutralino masses, of MSSM 
type, are modified to the order fi\/ fi K through the higgsino part, as is the 5th neutralino 
mass. The mass and the 55 wave-function are approximately given by 



r 



(0, 0, Sp, c p ) T 



(36) 




(37) 



and 




1 




(38) 



while doublet components are mixed in to first order, 



V*«-£?(0, 0, 8p , c p ) 



(39) 
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in parallel to the singlino components of the first doublet-type neutralinos. 

In summary, the gaugino/doublet higgsino dominated neutralinos follow the pattern 
of the MSSM quite narrowly. Increasing the value of \l k will increase the mass of the new 
singlino state (almost) linearly, causing the state to decouple and making the NMSSM 
very difficult to distinguish from the MSSM. 



3.3 The case with M x = M 2 in the limit of tan/? = 1 

When the two soft-breaking SU(2) and U(l) gaugino masses are equal, Mi = M2 = M 
and tan/3 = 1, cf. Ref.j^j, the electroweak gauge symmetry guarantees the existence of 
a physical neutral state which does not mix with the other states and which has a mass 
eigenvalue identical to the modulus M. Furthermore, the gaugino states do not mix with 
the singlino state S that couples only to the specific linear combination of the higgsino 
states H% = (#° + H G d )/^2. As a result, one gaugino state mixes only with one higgsino 
state while the other orthogonal higgsino state mixes with the singlino state, leading to a 
block-diagonal matrix composed of one scalar and two 2x2 matrices. 

This special structure can be made apparent by switching to the mixed basis 
{7, Z, H®, H®, S} from the original group basis {B, W 3 , H%, H®, S} by means of the trans- 
formation 



/ 7 \ 

Z 

Ha 

H b 

V s J 



A* 



W 3 
Hd 
H u 
\S ) 



( 



Cw Sw 

"% C\V 













1 

± 




\ 


-7 

* 

1/ 



w 3 

H u 
\S ) 



(40) 



In this new {7, Z, H®, H®, S} basis the mass matrix A4 5 takes the block-diagonal form 



Mb = AsMsAl 





M 





















M 


m z 















m z 


fi 



























V 


















(41) 



This mass matrix generates two two-state mixings between Z and H®, and between H® 
and S, respectively. The block-diagonal matrix can be diagonalized by the orthogonal 
matrix 



1 



V 5 



R 



g/h 



(42) 



R. 



h/s 



consisting of two 2x2 rotation matrices R g /h and R 



R 



g/h 



cos 9 g/h 
sin 9 g/h 



- sin 9 g/h 

COS 9 g/h 



R 



h/s 



COS 9 h/s 

sin 9 h/ s 



■ sin V h/s 
cos 9 h/s 



(43) 
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with the mixing angles determined by the relations 

2m z 



tan 0g/ h 
tan6 h / s 



M -/j- y/(M - p) 2 + 4m| 
2 fix 



+ fi+ \/ (jT k + /i) 2 + Aj? x 
The mass eigenvalues can be written completely in analytic form, 



mi = M 
1 



m 2 
m 5 



M + n+ \ (M — yu) 2 + 4m? 



fl K —fi—J (fl K + /i) 2 + 4yUl 



fi K -fi + J (ji K + /i) 2 + 4//1 



(44) 



(45) 



with the wave-functions of the neutralinos Xi [* = 1, . . . , 5] determined by the cos/sin of 
the mixing angles 6 g /h and 9h/ a - 
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Figure 3: (a) The neutralino masses \rrii\, mapped onto positive values, and (b) the 
tangent values of the mixing angles 9 g /h and 6h/ s as a function of the higgsino mass 
parameter \i for the parameter set: M = 200 GeV, \i K = 120 GeV, fi\ = 100 GeV. 

It is instructive to study the neutralino mass spectrum in this model for a set of fixed 
parameters M = 200 GeV and /i K = 120 GeV, ji\ = 100 GeV, by varying the higgsino 
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mass parameter /i. The branch character of the eigenvalues {23} and {45} is exemplified 
in Fig|3fa). The tan's of the corresponding mixing angles are displayed in Fig|3fb). With 
rising tan's we move from scenarios of no mixing, to maximal gaugino/doublet higgsino 
in {23} and doublet /singlet higgsino mixing in {45}, finally to gaugino/doublet higgsino 
flipping {23}, and doublet /singlet flipping {45}, while the gaugino {1} remains untouched. 



4 Neutralino Production and Decays 

In the MSSM the neutralino sector consists of two gauginos and two higgsinos. Typically 
the lightest supersymmetric particle (LSP), which is stable under the assumption of im- 
parity conservation, is the lightest state of the neutralino mass matrix. The LSP will 
appear as one of the final states of each sparticle decay and its non-observability is 
responsible for the well-known missing energy /momentum signature of supersymmetric 
particle production. 

The neutralino production and decay properties in the NMSSM with the additional 
singlino state depend crucially on the singlino mass with respect to the MSSM neutralino 
masses jS]. If the singlino is much heavier than the other states, it will be very rarely 
produced and so practically unobservable. On the contrary, if the singlino is lighter than 
the other states, a singlino-dominated state will be the LSP so that the other neutralino 
states will decay, possibly through cascades, into the singlino-dominated LSP. 

In this section, we present a qualitative description of the production of neutralinos, 
involving at least one singlino-dominated state, such as X5X5 an d X1X5 an d the subsequent 
decays of the neutralino Xi into leptons and light Higgs bosons. 



4.1 Singlino Production in e + e Annihilation 

The production processes 

e + e"^x°X? [i,j = l-5] (46) 

are generated by s-channel Z exchange, and t- and w-channel cl,r exchanges. After 
appropriate Fierz transformations of the selectron exchange amplitudes [with the electron 
mass neglected], the transition matrix element of the production process can be written 

as 

T{e + e- £ Q aP [v(e + ) 1 ,u(e-)] a [u{^v{x])] p (47) 

a,0 =L,R 

The transition amplitudes are built up by the sum of the products of chiral neutralino 
currents and chiral fermion currents. The four generalized bilinear charges Q a @ corre- 
spond to independent helicity amplitudes, describing the neutralino production processes 
for polarized electrons/positrons ^3]. They are defined by the fermion and neutralino 
currents and the propagators of the exchanged (s)particles as follows: 

QlL = + 2 ^2 (^3 ~~ Qf S w)^ij ~ DuLgLij, Q RL = 2~ Qf^ij + DtRdRij 

s w c w c w 
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Qlr = -- (4 - Q f s 2 w )Z:, + D tL g* Lij , Q RR = +^Q f Z* j - D uR g* Rij (48) 

with / = e~ in the production channel. The first term in each bilinear charge is generated 
by Z exchange and the second term by selectron exchange; D z , D tLR and D uLR denote the 
s-channel Z propagator and the t- and w-channel left /right-type selectron propagators 

s s 

Dz = n——. =r, D( t ,u)L,R = 77 — \ 2 — ( 49 ) 

s — m z z + irrizi- z [t, u) — m ~ 

J L ,R 

with s = (p e - +p e +) 2 , t = {p e - — p^o) 2 and u = (p e - —pyo) 2 representing the Mandelstam 
variables for neutralino pair production in e + e~ collisions. Finally, the matrices 2y, guj 
and g R ij are products of the neutralino diagonalization matrix elements Nf- 

g Lij = (ilN&cw + (Q f - I&Nfew) (l}N? 2 cw + (Q f - I^N^swJ /s 2 w c 2 w 
9m = Q)NlN^/c 2 w (50) 

The e + e~ annihilation cross sections follow from the squares of the bilinear charges, 
° [e + e- - = % ^A{£ £ { [1 - - n)f + A PS cos 2 6] Q 1 

1/2 "1 

+4fj Ji fj J jQ 2 + 2A P g Q 3 cos > dcos 9 (51) 

where <S^ is a statistical factor: 1 for % ^ j and 1/2 for i = j; fa = m^o / y/s, 6 is the polar 
angle of the produced neutrinos; and Aps = Aps(l, , faf) denotes the familiar 2-body 
phase space function Aps(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz. The quartic charges 
Qi (i — 1, 2, 3) are given by the bilinear charges as follows: 

Qi = \ [\Qrr? + \Qll\ 2 + \Qrl\ 2 + \Qlr\ 2 } 
Q2=\$e [QrrQ* rl + QllQIr] 

Qs = \ [\Qrr\ 2 + \Qll\ 2 - \Qrl\ 2 - \Qlr\ 2 } (52) 

An example for the production of xt i n association with another singlino-type (xs), 
or a gaugino-type (xT) or a higgsino-type (xl) neutralino is presented in FigJH for the 
parameter set IP [FigJT] with mg R = 200 GeV and mg L = 250 GeV. [Of course the {55} 
final state is unobservable without additional ISR 7 emission.] The increase of the cross 
sections with increasing doublet-singlet gaugino/higgsino mixing parameterized by fix is 
obvious. The gaugino character of Xi is responsible for the dominant size of the {51} 
cross-section. 

With the anticipated integrated luminosity f £ = 1 ab _1 , sufficiently large event rates 
of order 10 3 are predicted if n\ is not too small. 
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% [GeV] 

Figure 4: The production cross sections of neutralino pairs, {51} (dashed), {55} 
(thin- s olid) , {53} (dotted) and {11} (thick-solid), in e + e~ collisions with the center- 
of-mass energy y/s = 500 GeV as a function of fi\ for the parameter set P [FigU^ with 
m iR = 200 GeV and m e ~ L = 250 GeV. 



4.2 Decays to a Singlino, with no Higgs bosons 

(i) If kinematically allowed, two-body decays of neutralinos to the electroweak gauge 
bosons Z are among the dominant channels. The widths of decays Xi ~^ X*jZ are given 
by 



n 2 \ 1/2 

1677771^0 

Xi 



\z 



l 3 



(m| -m|o) 2 

— n — hm|o+m? -2m| 

mi x * K i 



+ 6m^oms-o^te(Z^) > (53) 



where Aps = Apsfl, m 2 ~ /m1 , m 2 z /m 2 ~ ), with Zij defined in Eq. (J5UJl . The widths of the 

Xi ' v 1 



chargino 2-body decays into a neutralino and a W boson, \t ~ > Xj W^S rea d correspond- 
ingly 



g A PS 
167rm-± 

Xi 



+ |w 



(m.± — m?q) 

Xl 2 Xj + m'± + m 2 - 2m 2 

Xi Xj 



m 



w 



-6 m-±m^o afte(W W j W^) } 



(54) 



where Aps = \ps(l, m 2 -q/m 2 ± ,m 2 v /m 2 ± ) and the bilinear charges Wl,r are defined as 



Xi 



U* Rll N% - ±=U* R ^l (55) 
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The unitary matrices Ul and Ur diagonalize the chargino mass matrix as 
U R M c Ul = diag|m-±,m-±|, cf. Ref.[H| for details. 

If 2-body decay channels are closed kinematically, the 3-body neutralino decays, 
Xi ~^ xlffj are generated by s-channel (virtual) Z exchange, and t- and w-channel 
sfermion exchanges. Neglecting fermion masses, the transition matrix element, cf. 
Ref.|15j. is determined by the bilinear charges Q' a/3 which are related to the bilinear 
charges Q a p introduced for the production, by crossing symmetry as 

4aP = Qlp ( 56 ) 

with the transformed Mandelstam variables, s = (pf + Pf) 2 , t = (p^o + pj) 2 and 
u = (p^o +pj) 2 for the decays. [Neutralino decays to charginos and W bosons can be 
described in the same way after obvious redefinitions of the bilinear charges.] Decay 
widths and distributions depend on the quartic charges Q[, Q' 2 and Q' 3 defined analo- 
gously to Eq. (J52J) . 

(ii) At the LHC, cascade sfermion decays, / — > /%°, are of great experimental interest. 
The width of the sfermion 2-body decay into a fermion and a neutralino follows from 

r|/ - fi ' ] = w^H~ f H " m h - '"/) < 57 > 

where the 2-phase space function Aps = Aps(l, m 2 o/m 2 f , m 2 f /m 2 f ) with / = f L or f R ; the 
couplings are expressed in terms of the neutralino mixing matrix A^ 5 as 

g hi = V2 [l(Nf 2 + (Q f - 4)Nlt w ] and g f ~ Ri = 2Q f s w N?* (58) 

in obvious notation. 

The reverse decays, neutralino [chargino] decays to sfermions plus fermions, Xi ~^ ff 
etc, are given by the corresponding partial widths, 

rR - - //] = £&o ls /« |2 H + m < - "*/) (59 > 

Xi 

with the same couplings as before and Aps = Aps(l, m 2 Jmio, m 2 /miq). [Analogous ex- 
pressions hold for chargino decays.] 

Examples of these partial decay widths are shown in Fig|S] [with the parameter set P 
as in FigJTj. For an illustrative purpose, the mass of the R-sleptons Ir is assumed to be 
mj R = 200 GeV > m^o for the 3-body neutralino decays and to be m- [R = 130 GeV < m^o 
for the 2-body slepton decays, respectively. § The masses of the squarks are assumed to 
be rriq L = 250 GeV and rriq R = 200 GeV. For small mixing fi\ the lifetimes of the second 
lightest neutralino xi an d the R-sleptons Ir can be quite large, giving rise potentially 

§In either case Xi or Ir is the next-to-lightest SUSY particle NLSP with just one decay channel open 
to the lightest SUSY particle LSP= xi- 
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to macroscopic flight paths 9 . However, cosmological bounds on must be analyzed 
before any (realistic) experimental conclusions can be drawn. The kink in the Xi lifetime 
and flight distance in the upper right panel of Figgis caused by accidental cancellations 
between sfermion and Z exchange diagrams in the decays Xi ~^ X5 QQ an d X5 y v] these 
accidental cancellations do not occur [to any significant degree] in the decay Xi ~~ y xl 




\i, [GeV] 




\i x [GeV] 




\i x [GeV] 



Figure 5: The widths, lifetimes and flight distances [broken lines] of the decays Xi ~ > 
%5 and Ir — > as a function of fi\ for the parameter set FfFigU^. The mass of 
the right-handed slepton is taken to be m~ lR = 200 GeV > m^o for the 3-body neutralino 
decays (upper panels) and to be m~ lR = 130 GeV < m^o for the 2-body slepton decays (lower 
panels). The masses of the squarks are assumed to be rriq L = 250 GeV and rriq R = 200 
GeV. Right: flight distances for a/s = 500 GeV are shown by broken lines. The kink in the 
Xi lifetime and flight distance in the upper right panel is caused by accidental cancellations 
between sfermion and Z exchange diagrams. [The value of the lower bound, expected from 
cosmological arguments on //> is presently not yet known.] 
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4.3 Decays to a Singlino, involving Higgs bosons 

Decays involving Higgs bosons can be quite different for different Higgs boson mass spec- 
tra. Following the procedure outlined in Ref. 0] we decompose the neutral Higgs states 
into real and imaginary parts as follows: 

H° d = -^=[v d -S 1 s f3 + S 2 c f3 + tP 1 s p }, (60) 
H° u = -^=[v u + S l c p + S 2 sp + iP 1 cp] } (61) 

S = -j=[v a + S 3 + iP 2 ] (62) 

where the Goldstone states are removed by using the unitary gauge. We then further 
rotate these states onto the mass eigenstates, Hi (i = 1 — 3) and A{ (i = 1, 2) labeled in 
order of ascending mass, by using the orthogonal rotation matrices^ H and A : 

Hi = S 3 0% A, = pp* (63) 

The resulting mass spectrum, composed of three scalars, two pseudocalars, and two 
charged Higgs bosons, is shown in FigE] as a function of For the purposes of ex- 
ample, we have chosen the mass parameter Ma (defined to be the heavy pseudoscalar 
mass in the MSSM limit) to be 2/i/ sin 2(3 ~ 567 GeV, setting the scale of the heavy 
Higgs bosons. The lighter Higgs bosons consist of two scalars and one pseudoscalar. The 
lightest scalar and pseudoscalar in our example are predominantly singlet states, with 
masses set by the scale of \i K . 

Generally, the width of a 2-body neutralino or chargino Xi decay to a neutralino or 
chargino Xj an d a Higgs boson <pk (Hk or Ak) is given by 

m~m = is|r { K + «• - <) <w + W) 

Xi 

+2 W m^m*. (Cg fc Cg; + CgfCgO } (64) 

where Aps = Aps(l, m| /m|., m^/m-J and the left/right couplings C\j k R must be speci- 
fied in each individual case; rj^ = 1 for <fi = H^, H ± , and —1 for <ft = A*. 

(i) For the decay of a neutralino Xi to a neutralino x°j an d a scalar Higgs boson Hk, 
Xi — > x'jHk, the couplings are given by, 



Cg k (Xi ^ $ H k) = l[g(N! 2 -Nlt w )( Nf 3 sp + N^) + V2\(N? 3 cp-N? 4 sp) N; 



2 
1 

+ 2 

+ ^[XN^-KN^Og + (65) 



^Note that the definitions of these mixings matrices differ slightly from those in Ref.0], where the 
scalar rotation matrix is defined via O t — H and the pseudoscalar rotation is defined by a rotation 
through an angle 6 a- 
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Figure 6: The Higgs boson mass spectrum as a function of fix for the parameter set 
P/FzgOF an d maximal mixing. For the purposes of example, the Higgs mass parameter 
Ma is set to 2/x/ sin 2f3. Heavy scalar, pseudoscalar and charged states are nearly mass 
degenerate: M Hs ~ M A2 ~ M H ± ~ M A . 



While the first term in each of the two square brackets in Eq. (Jo5|) are reminiscent of the 
MSSM couplings XiX'jh and x°iX°jH respectively, the other terms are genuinely new in 
origin, arising from the extra interaction terms in the NMSSM superpotential. 

The widths for the kinematically allowed decays x\ — > Xi 5-^1,2 are shown in Fig0 
(left) as a function of fix- For fix = the xt state is decoupled from the other neutralinos; 
as fix is switched on, the coupling, and therefore the decay widths, increase. The decay 
widths for xi — ► X5H1 an d xi ~~ > X5H2 are comparable, within an order of magnitude, due 
to the large xi mass and the near mass degeneracy of Hi and H 2 . With partial widths of 
order GeV, these decay modes are in the observable range of branching ratios. 

(ii) Similarly, a 2-body neutralino decay to a neutralino and a pseudoscalar Higgs 
boson, x° — > x^jAk, follows Eq. ljMjl with the left/right couplings given by 



Cijkix'i 



xjA k ) 



+ 



g{Nl-Nlt w ){-N%s p + N^c p ) + V2\ (N? 3 c p +N? 4 s p ) Nf, 
-^=[XNlN^-KN! 5 N^}Oi + 



(67) 
(68) 



Again, only the first term in the square brackets is similar to the MSSM coupling XiXj-^-- 

The widths for the kinematically allowed decays xi ~^ Xi 5^1 are shown in Figd (right) 
as a function of fix for our chosen example scenario. In comparison with the scalar case, 
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Figure 7: The decay widths for x® — » 5 Hi 2 (left) and x° - ► Xi 5^1 (right) for the 
parameter set ¥[Fig\^j and maximal mixing. For the purposes of example, the Higgs mass 
parameter Ma is set to 2\xj sin 2f3. 

many of the decays are kinematically disallowed, only leaving the decays of the heaviest 
two neutralinos to xl an d the lightest pseudoscalar (Ai). Note that pseudoscalar decays 
are strongly suppressed compared with the scalar modes and may not be observed easily. 

(iii) For completeness, we describe the decays of charginos to a neutralino and charged 
Higgs boson xf ~ > Xy,H ± (i = 1,2; j = 1 — 5). These follow a similar pattern, now with 



the last index of the coupling removed: 



-^(N? 2 *+N?*t w )U* Lj2 



N 5 *TT* 
iv i4 U Ljl 

1 



-\s p Nl*Ul j2 (69) 
-\ Cf3 N? 5 U* Rj2 (70) 



However, the large mass of the charged Higgs boson means that these 2-body decays are 
kinematically disallowed for our specific parameter choice. 

(iv) It is also possible for Higgs bosons themselves to decay into the singlino- 
dominated state, via the decays Hi — > xix^ M — > x^>X°j an d H ± — > xixt^ if kinematically 
allowed. Clearly this is only possible for the heavier Higgs states; the lightest Higgs boson 
is never heavy enough to decay in this way. The general form of the width for these decays 
4>i ~^ XjXk (<Pi = Hi, Ai, H^), is given by the crossing of Eq. lfMJ) : 

,V2 



-2 V4> m^ k ((%..(%: + C^ k Cf jk ) } (71) 

where Aps = Aps(l, m|./m^, rri^./m^.) and Sjk = 1 or 1/2 is the usual statistical factor. 
Again, rj^ = 1 for <p = Hk, H ± , and — 1 for cf) = A k . The couplings C\j k R are related to 
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their neutralino decay counterparts in the obvious way: 

C^H^xUl) = C^ix , ^ X° k Hi) (72) 

C^iA^m) = ^(X^XIA) (73) 

<% ,R {H± - xUf) = C^ixf^x ^) (74) 

Some of these decays widths are plotted in FigJHJ Note that a significant fraction of the 
Higgs boson H% and A 2 decays go into the invisible channel X5X5 only if the partial decay 
width exceeds the range of ~ 1/10 GeV. The upper left panel, showing the partial width 
for the decay H 3 — > X4X5 [i — 1; 5], has been allowed to extend down to widths of order 
1CT 5 GeV to show the switching off of the X1X5H3 coupling at ~ 58 GeV. This is caused 
by destructive interference between the different constituent fields in both the Higgs and 
the neutralinos, and is directly analogous to the cancellations seen in Ref. jl] . 




Figure 8: The decay widths for H 3 — > Xi 5X5 ( u PP er left), A 2 Xi ; 5X.5 (upper right) 
and — > xtx.i 5 (lower) for the parameter set F[Fig^ and maximal mixing. For the 
purposes of example, the Higgs mass parameter Ma is set to 2fi/ sin 2f3. 
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5 Summary and Conclusions 



In this study, we have investigated the neutralino sector of the NMSSM, suggested by 
many GUT and superstring models. Moreover, this model attempts to explain the \x- 
problem of the MSSM by introducing a new iso-singlet Higgs superfield, S, the scalar 
component of which acquires a non-zero vacuum expectation value. 

We have given expressions for the new 5x5 neutralino mass matrices and mixing 
matrices and we have presented, besides the numerical analyses, approximate analytical 
solutions for the neutralino masses and mixings which provide a nice insight into the 
structure of the spectrum and the mass hierarchies in case of small couplings between the 
MSSM and the new iso-singlet. 

The renormalization group flow of the parameters A and k from the GUT scale down to 
the electroweak scale gives rise to strong upper bounds on their values at the electroweak 
scale, where small k is favored. The qualitative features of the neutralino masses are 
dependent on how strongly the PQ symmetry of the model is broken by non-zero k 
values; this is quite accurately described by the approximate analytical solutions. 

If the PQ symmetry is slightly broken for small k, the qualitative pattern for the 
particle spectrum remains intact, except that the lightest singlino- dominant neutralino 
acquires a mass of the order of the electroweak scale. Thus the model contains four 
MSSM-type heavy gaugino/higgsino dominant states and one light singlino dominant 
state. Since the couplings to the Z boson can be very much reduced, the NMSSM with a 
slightly broken PQ symmetry constitutes a valid scenario. 

In contrast, a strongly broken PQ symmetry, though disfavored by the flow of the 
couplings from the GUT scale down to the electroweak scale, could provide an extra 
moderately heavy neutralino state, which is only weakly coupled to the Z and (s)fermions. 
Such decoupled scenarios would be more difficult to distinguish the NMSSM from the 
MSSM. 
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Appendix: The small— mixing approximation 

The 5x5 neutralino mass matrix of Eq.(j3J in general cannot be diagonalized analytically to 
derive the physical neutralino masses. However, simple analytical expressions for masses 
and mixing parameters can be found by making use of approximations for small doublet- 
singlet mixing which is theoretically very well motivated. 

To construct this approximate solution in the neutralino sector, we treat the doublet- 
singlet mixing parameter fi\ <C Mj,/i, together with the Z-boson mass mz, as small 
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parameters of generic size e <C 1 in units of the typical SUSY masses. Then, as long 
as these SUSY masses are not as small as fix, we observe a hierarchical structure in the 
neutralino mass matrix of the form: 

« = ( d n) ( 75 ) 



x X B , 

where A is a 4 x 4 matrix incorporating elements of the order of the large SUSY scale, B 
is a scalar and X is a 4-component vector of order e. 

Performing an auxiliary orthogonal transformation O defined by the matrix" , 

0=\ 2 , 76 

with the mixing column vector Q = [A — B]~ l X, the mass matrix takes block diagonal 
form, accurate to order e 2 : 

rp ( A + A A \ 

0H0 T = J 77 

V B + A B J K J 

where 

A A =1{[A-S]- 1 ,X^} 

A B = -X T [A- B]~ 1 X (78) 

Both A's are of order e 2 . If A is diagonal, only the diagonal elements of need be 
kept as re-diagonalization would change the mass matrix (|77|) and the orthogonal matrix 
(|76|) only beyond the order considered in the systematic expansion. We note that the 
correction terms satisfy the simple sum rule TrA^ + A# = 0. 

If B is also as small as the elements of the low vector X, the mixing column vector 
Q = A~ X X and the correction terms A^ and A# are further simplified to be 

A A = - {A~\ XX T } and A B = -X T A~ X X (79) 

On the contrary, if B is much larger than the other parameters, the mixing column vector 
Q = —X/B and the correction terms take the following simple form 

A A = -XX T /B, A B = X T X/B (80) 

Both these approximations have been used in the derivation of all the mass and mixing 
formulae discussed earlier in the report. 

The diagonalization of the mass matrix A4§, 

(8l) 



"Note that by standard notation QQ T is a 4 x 4 matrix with the elements QiQj while fl T fl is a scalar 
with the value YJ Qf. 
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with M. being the 4x4 MSSM mass sub-matrix, makes use of the block-diagonalization 
method in the following way: 

(1) In the first step M. is diagonalized by applying the well-elaborated MSSM procedure 



M u = VMV 1 (82) 
generating the eigenvalues M. D = diag [mi, . . . , m 4 ]. 

(2) The ensuing 5x5 matrix can subsequently be block-diagonalized as worked out above 
by applying the orthogonal transformation in Eq. (|76j) with 

n = VT : r = V T {M D - m^' 1 V X = (M-m^X (83) 

Note that T is of order e - quantum mechanically enhanced however if mass differences 
{"ffn — are, accidentally, small. 

(3) The block-diagonalization affects the upper left diagonalized 4x4 submatrix Ai D 
only beyond second order and likewise the orthogonal matrix V 5 beyond the second and 
first order considered, respectively, for on- and off-diagonal elements. As a result, we 
obtain the final diagonal form of the mass matrix as 

Mg « V 5 M 5 V 5T (84) 

with 

/wKvneviT CVT) \(v o \ 

V -(VTf 1 - ±(W) T (VT) ) \ l lxl j 

in obvious notation. While the right-most part solves the MSSM diagonalization, the 
left-most part diagonalizes the NMSSM under the assumption of small doublet-singlet 
mixing. 
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